function dx=Plant(t,x,flag,para)
dx=0;

Zeta=para(1);
dZeta=para(2);

tf=2;delta=0.00001;
kg=dZeta/(1-Zeta+delta);

dx(1)=-kg*x(1);  
